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It has become clear in the past several years that the cuprates show many unusual properties, 
both in the normal and superconducting states, especially in the underdoped region. In particu- 
lar, gap-like behavior is observed in magnetic properties, c-axis conductivity and photoemission, 
whereas in-plane transport properties are only slightly affected by the pseudogap. I shall argue that 
these experimental evidences must be viewed in the context of the physics of a doped Mott insu- 
lator and that they support the notion of spin charge separation. I shall review recent theoretical 
developments, concentrating on studies based on the t-J model. I shall describe a model based on 
quasiparticle excitations, which predicts the doping dependence of T c and anomalous energy-gap- 
to-T c ratios. Finally, I shall outline how the model may be derived from a microscopic formulation 
of the t-J model. After a brief review of the U(l) formulation, I shall explain some of the difficulties 
encountered there, and how a new SU(2) formulation can resolve some of the difficulties. 



I. INTRODUCTION 

It has become clear in the past seveal years that the 
cuprates show many highly unusual properties both in 
the normal and superconducting (SC) states. These un- 
usual features are related to the fact that the cuprates 
are doped Mott insulators. It is then not surprising that 
the unusual behaviors are most striking in the under- 
doped region, when the concentration of doped holes, x 
is small. In the normal state a pseudogap is obseved in a 
temperature range considerably above the SC transition 
temperature T c . The gap is seen in NMR relaxation rate 
1/Ti, Knight shift |l|] and specific heat. Q It is also seen 
in c-axis conductivity ||] and in photoemission experi- 
ments ||Jj| which reveal that the pseudogap is roughly 
of the same size and k depdendence as the d-wave SC 
gap. Furthermore, the gap size is essentially indepen- 
dent of x and even increases slightly when T c is reduced 
with decreasing x. This observation is also supported by 
tunneling data. Q On the other hand, the in-plane trans- 
port properties are only slightly affected by the pseudo- 
gap. The resistivity shows a small decrease which may 
be interpreted as a decrease in the scattering rate. 0] 
More importantly, the spectral weight of the Drude part 
of <j{lo) is proportional to x Q and there is no evidence 
that it is strongly reduced by the presence of the pseu- 
dogap. We believe this is strong experimental ev- 
idence supporting the notion of spin-charge separation 
JTl| in these materials. It was pointed out by P.W. An- 
derson early on that the Neel state is not the best 
way to accommodate the competition between the hole 
kinetic energy and the spin exchange energy. He envi- 
sioned another possibility, i.e., the spins form a liquid of 
singlets, which he termed the resonating valence bond 
(RVB) state. The reason is that the energy to form a 
singlet —JS(S + 1) is particularly favorable for S = \. 
The holes can more freely among the liquid of singlets 
and are responsibile for the charge transport. This no- 
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FIG. 1. A cartoon representation of the RVB liquid or 
singlets. Solid bond represents a spin singlet configuration 
and circle represents a vacancy. In Fig. (lb) an electron is re- 
moved from the plane in photoemission or c-axis conductivity 
experiment. This necessitates the breaking of a singlet. 



tion of spin charge separation naturally accounts for all 
the qualitative features of the spin gap state noted above. 
The spins form RVB singlets so that it costs energy (spin 
gap) to make triplet excitations. However the in-plane 
conductivity is carried by x holes, which remain gapless. 
In c-axis conductivity and photoemission, a physical elec- 
tron is removed from the plane, which carries both spin 
and charge. It then follows that the spin gap should ap- 
pear in these experiments. This picture is illustrated in 
Fig. 1. 

We note that an alternative model which exhibits the 
above phenomenology is the model of preformed pairs 
above T c . There are two versions of this class of model; 
the first suggests that strong phase fluctuation Jl^] de- 
stroys long range order over a large temperature range 
above T c , and the second assumes that we are in the short 
coherence length limit of the pairing state , so that es- 
sentially pair "molecules" are first formed and then Bose 
condensed. JUIyl The phase fluctuation model predicts 
that the pairing amplitude is responsible for the pseudo- 
gap and would seem to predict that other manifestations 
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of superconductivity such as conductivity and diamag- 
netism fluctuations should be observable, particularly at 
short distance and short time scale. A recent high fre- 
quency conductivity experiment in underdoped BISCO 
p6[ shows that while Berenzinskii-Kosterlitz-Thouless 
(BKT) type fluctuations are observed near and above 
T c , the short distance (bare) superfluid density extracted 
from these measurements vanishes above 100 K, much 
below the temperature range associated with the pseu- 
dogap. These data are difficult to understand within the 
phase fluctuation model. Similarly, in the short coher- 
ence length model, charge transport is by charge 2e pairs 
in the pseudogap state and it is difficult to understand 
the insensitivity of the transport properties to the ap- 
pearance of the pseudogap with underdoping. Further- 
more, it is not at all clear that the coherence length is 
short in the underdoped limit. In section III we shall in 
fact argue that the coherence length increases with un- 
derdoping, and that one is not in the short coherence 
length regime. In any event, in both these models a su- 
perconducting state with a large energy gap is postulated 
to exist, without any indication of the origin and the en- 
ergy scale of the gap. The RVB picture is fundamentally 
different from these preformed pair pictures in that spin- 
charge separation plays a crucial role. The pseudogap is a 
spin gap with an energy scale set by J, which becomes the 
superconducting gap with the onset of coherence in the 
charge degrees of freedom. The superconducting state 
is characterized by spin-charge recombination, forming 
superconducting quasiparticles which are quite conven- 
tional in the BCS sense. 

We model the cuprate with the t-J model, which we 
believe contains the essential physics of the doped Mott 
insulator. The t-J Hamiltonian is 
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is subject to the constraint that double occupancy of a 
site by two electrons of opposite spins is not allowed. 
Here Si = c\ a cr a pCij3 and rn = J2 a c la c ia- Tne t-J model 
is the strong coupling limit of the Hubbard model and the 
difficulty of its solution lies in enforcing the no double 
occupancy constraint. For the cuprates the parameters 
are known to be J ~ 0.13 meV and t/J w 3. When 
holes are doped into the insulator, there is a gain in ki- 
netic energy per hole proportional to t due to hopping. 
However, the spin correlation is destroyed, costing an en- 
ergy of approximately J per site. Thus we can consider 
the doping problem as a competition between the en- 
ergy xt (kinetic energy per site) and J. When xt « J, 
the AF state with its doubled unit cell is retained and 
the holes form small pockets around the top of the sin- 
gle hole dispersion, which is known to be at (vr/2, ±7r/2) 
from photoemission [jl7j (see Fig. 2). This problem be- 
longs to the same class as the doping of a band insulator 
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FIG. 2. The evolution of the locus of low lying single parti- 
cle excitation with doping, (a) low doping AF state: Brillouin 
zone is doubled with small hole pockets of area x. (b) un- 
derdoped: Fermi surface "segments." (c) overdoped: Fermi 
surface with area 1 — x, satisfying Luttinger theorem. 



(or semiconductor). The only difference is that the co- 
herent part of the band has a reduced spectral weight 
of J/t and a bandwidth of order J. This can be under- 
stood in terms of a spin-polaron picture, i.e., the hopping 
hole is surrounded by a cloud of spin excitations. [H] On 
the other hand, it xt » J, the spin correlation becomes 
unimportant, AF order is destroyed, and the holes should 
form a metallic state, describable by Fermi liquid theory. 
The Luttinger theorem then dictates that the area of the 
Fermi surface is given by 1 shown in Fig. 2c. The 

important point is the electrons which form the local mo- 
ments in the Mott insulator are now mobile and should 
be counted as part of the Fermi sea. The change in Fermi 
surface area from x to I — x between the low and high 
doping limit is a special feature of the doping of a Mott 
insulator, associated with the liberation of the local mo- 
ments. The question then arises: how does the system 
evolve between these two limits? The intermediate state 
is apparently the spin gap state, with gaps in the one 
electron spectrum in the vicinity of (0, n) and segments 
of the Fermi surface near (7r/2,7r/2) (see Fig. 2b). As 
doping is increased, these segments grow in length and 
eventually join to form the Luttinger Fermi surface. This 
intermediate state is clearly not a Fermi liquid because 
in band theory, gapping of parts of the Fermi surface is 
not permitted without symmetry breaking. However, the 
breakdown of Fermi liquid theory is not a sharply posed 
issue at finite temperatures. The existence of the super- 
conducting ground state at intermediate doping means 
that this question cannot be investigated experimentally 
at present. It is worth noting that the transition region 
occurs near x = 0.2, when xt and J are comparable. 

We emphasize once again that an important aspect 
of the doping of the Mott insulator is that the result- 
ing metal must remember that x holes are responsible 
for the electrical conductivity. If the AC conductivity 
cr(u>) is characterized by a Drude-like component at low 
frequency, we may characterize the conductivity by the 
scattering rate 1/r and the spectral weight (n/TO) e ff e ctive- 
For underdoped samples, this spectral weight is propor- 



2 



tional to x. || This is very natural in that the weight 
must vanish when x — > 0. It follows that a superconduc- 
tor that forms out of the underdoped metal must have a 
superfluid density p s given by this spectral weight (in the 
clean limit), so that p s is proportional to x. This simple 
observation will play a prominent role in our subsequent 
discussion. On the other hand, when xt > J we have 
a Fermi liquid state with electron density 1 — x. The 
question then arises as to how (n/m) e ff e ctivc ~ % can be 
accommodated within Fermi liquid theory. Within Fermi 
liquid theory we can write 
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where to* is the effective mass and F\g is a Landau pa- 
rameter. It describes the deviation of the current carried 
by the quasiparticle from — evt due to dragging of other 
quasiparticles 
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where a=(l + Fig/2). |Q Here we have made the sim- 
plifying assumption that only the 1=1 Landau parame- 
ter (in 2d) is important and the correction is independent 
of k. From Eq. (1) we see that there are two ways to ob- 
tain (n/m) e ff = x. The first is to generate a heavy mass 
so that to* — 1/x. This is in fact the case for the system 
Lai^^Sr^TiOa which is a Mott insulator for x — with a 
Neel temperature of T N = 150 K. With doping x > 0.02, 
a metallic state is formed with the spin susceptibility \ 
and the specific heat coefficient 7 both scaling as 1/x 
and a Wilson ratio of order unity. p(| The Hall coeffi- 
cient i?jv ~ 1/(1 — x) as expected for a Fermi surface 
dictated by Luttinger theorem. This is clearly a realiza- 
tion of the Fermi liquid state expected for xt > J. Unlike 
the cuprates, LaTiC>3 is a three dimensional system, so 
that the exchange constant J can be deduced from the 
ordering temperature. Thus the ratio J/t is very small 
and we believe this is the reason why the Fermi liquid 
state persists to low doping. For x < 0.02, disorder ef- 
fects become important and we are not able to explore 
the xt < J limit in this system. 

A second route to achieve a spectral weight of x is for 
(l + ^A 5 -) » x. It turns out this is the route followed by 
the mean field slave boson theory described below. 
We shall see that in the underdoped region this route 
is not followed in the cuprate system. We have strong 
evidence that the factor a in Eq. (2) is not proportional 
to x and is in fact of order unity. 



II. MICROSCOPIC MODEL AND MEAN FIELD 
THEORY 

The physics of spin charge separation appears natu- 
rally in a class of theory which starts with the t-J model 



and enforces the constraint of no double occupation by 
decomposing the electron into a fermion and a boson, 
c\ a = ]\ a bi- The fermion f ia carries spin index and 
the boson 6j keeps track of the charge degrees of free- 
dom. The constraint is replaced by the requirement that 
fiafiiy + b\bi = 1 which can be enforced by introducing 
a Lagrangian multiplier so that field theoretic methods 
may be applied. This decomposition (called the slave 
boson method) is not unique and one could just as well 
associate the spin with the boson (the Schwinger boson 
theory |2^] ) . If the theories are solved exactly they should 
give identical results. However, different factorization 
leads naturally to different approximation schemes. Our 
strategy is to explore the different schemes to see which 
correspond most closely with experiment. In particu- 
lar, while the Schwinger boson method gives an excellent 
description of the antiferromagnetic state at half filling, 
|p2| it does not produce a large Fermi surface for large 
doping. Since we are mainly interested in the regime of 
intermediate doping, the slave boson is a more promising 
starting point. 

The exchange term can be written in terms of the 
fermions only j23[ 
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which invites the following mean field decoupling 
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and 



Ay - (fitfji - filfjf) 



(4) 



These parameters describe the formation of a singlet on 
the bond ij. The mean field phase diagram |24|,|2^] is 
shown schematically in Fig. 3. As the temperature is 
lowered, Xij 7^ 0, so that the fermions now acquire an en- 
ergy band and a Fermi surface. At a lower temperature, 
the fermions form a pairing state with <i-wave symmetry. 
The bosons become essentially Bose condensed (with ex- 
ponentially large correlation length with decreasing T) 
below a cross-over temperature Tgg = 2irxt. Below Tg^ 
the boson field can be treated as a c-number. In the 
over doped region this gives rise to a Fermi liquid phase, 
similar to the theory of heavy fermion systems. In the 
intermediate doping range, the simultaneous presence of 
Ajj and < b > gives rise to a pairing order parameter for 
physical electrons (cijCjj,) which is of d-wave symmetry. 
Above T [ gl spin charge separation occurs at the mean 
field level. In the pairing state a <i-wave type gap occurs 
in the spin excitation spectrum, but not in the charge 
excitation, and it is natural to identify this as the spin 
gap phase. Finally, the region IV in Fig. 3 is a non 
Fermi liquid state which may be referred to as a "strange 
metal." 
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FIG. 3. Schematic mean-field phase diagram of the t-J 
model. Below the dashed line the uniform RVB order param- 
eter \ij is nonzero. The mean-field pairing line below which 
Aij 7^ (dotted) and the Bose-condensation line (solid) divide 
the phase diagram into four regions. Region I is a Fermi-liquid 
phase, region II is the spin-gap phase, region III is the super- 
conducting phase, and region IV is the strange metal phase. 

We can go beyond mean field to include fluctuations 
about the mean field solution. The most important fluc- 
tuations are the phase fluctuations of the order parame- 
ter Xij = \Xij . Particles hopping around a plaquette 
acquire a phase related to just like electrons in the 
presence of a magnetic flux. These low lying excitations 
are U(l) gauge fields. ]2(| We shall refer to this theory 
as the U{\) formulation. When coupled to the fermions 
and bosons they enforce the constraint locally, not just 
on average as in mean field theory. 

III. PHENOMENOLOGICAL DESCRIPTION OF 
THE SUPERCONDUCTING STATE 

Before continuing with the microscopic theory, we di- 
gress to review a phenomenological description of the su- 
perconducting state. J27]] The idea is to start at low tem- 
perature where the nature of the elementary excitations 
is well understood, and calculate the reduction of the su- 
perfluid density with increasing temperature. This can 
be done by making two assumptions: A) the superfluid 
density is given by x, and B) the quasiparticle (qp) dis- 
persion in presence of an external electromagnetic gauge 
potential has a BCS form: 

i4 sc) (k) = ±Ve 2 (k) + A2(k) + I ■ j(k) (5) 

for k near the nodes, where j(k) is given by Eq. (2). 
Note that v> = d^e is the "normal state" Fermi velocity 
and that the vector potential A couples only through 
the "normal state dispersion" e(k) and has nothing to 



do with the SC gap A(k). The physical reason for this 
is that this quasiparticle is a superposition of an electron 
with momentum k and a hole with momentum — k, and 
both these objects carry the same charge current j(k). 
Mathematically Eq. (5) is easily derived by noting that 
A enters only in the diagonal elements of the BCS matrix 
in the form e(k + A) — p, and — e(— k + A) — fi, which is 
diagonalized to give Eq. (5) to linear order in A. 

With these assumptions we can calculate how the su- 
perfluid density is reduced by thermal excitation of quasi- 
particles. We found that 

p. x _2^1n2^ T 

where V2 = Aoa/v2 is the velocity of the quasiparticle 
in the direction of the maximum gap Ao, i.e., in the di- 
rection from the node towards (0, tt). The ratio vpjv2 
thus measures the anisotropy of the massless Dirac cone 
which characterizes the d-wave qp spectrum. 

Based on numerical calculations and theoretical con- 
siderations, ]28| , |29t we expect the mass m in Eq. (6) to 
correspond to a tight-binding hopping integral of order 
t so that m~ x ~ t. Experimentally it is found that m 
is about twice the electron mass which happens to 
correspsond to a hopping integral of J ~ 0.13 cV. We 
believe that the theoretical expression for the hopping 
integral is i/3 which happens to equal J in our case. On 
the other hand, the Fermi velocity vf is proportional to 
the coherent bandwidth, which is given by J. To keep 
our expression general, we keep track of the distinction 
between i/3 and J, even though numerically they are 
equal. 

We see that for small x, the quasiparticle excitation is 
an effective way of destroying the superconducting state 
by deriving p s to zero. By extrapolating Eq. (6) to p s — 
0, we can estimate T c as being of order 

T c «xA (t/3J) . (7) 

Note that this is a violation of the BCS ratio 2A /kT c = 
constant. Presumably the real transition is driven by 
critical fluctuations including phase fluctuations and vor- 
tex unbinding in the 2d limit, but the underlying (bare) 
superfluid density should be driven to zero by quasipar- 
ticle excitations in the way we indicate. There is experi- 
mental support for this point of view from high frequency 
conductivity measurements. | fl6| If we further assume 
that Ao is independent of x for underdoped cuprates, 
we see that T c is proportional to x (or more precisely 
to p s (T = 0)/m), thus providing an explanation of Ue- 
mura's plot. Jl4| 

It is worth noting that in the slave boson mean field 
theory A is proportional to J, since that is the only 
energy scale relevant to the formation of spin singlets. 
Then Eq. (7) predicts that T c is proportional to xt. Apart 
from numerical coefficients, this has the same functional 
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dependence as the Bose condensation temperature T^X, 
as well as the transition temperature based on pairing of 
bosons to be discussed later in the SU(2) formulation. 

Another important implication is that superconductiv- 
ity is destroyed when only a small fraction of the quasi- 
particles (with energy < xAq) are thermally excited. 
Thus the gap near (0, 7r) must remain intact in the normal 
state, leaving a strip of thermal excitations which extend 
a distance proportional to x from the nodal points. This 
is qualitatively in agreement with the photoemission ex- 
periment. Of course our phenomenological picture does 
not provide a description of the normal state. It simply 
states that the normal state gap is an inescapable conse- 
quence of a finite Ao and a vanishingly small superfluid 
density as x — > 0. 

The fact that dp s /dT is independent of x and that both 
p a and T c are proportional to x means that a scaled plot 
of Ps{T) I p s {0) vs T/T c should be independent of x for 
small T/T c . In fact, such a scaled plot for YBCO6.95 and 
YBCO6.60 shows a remarkable universality over the entire 
temperature range. Q We can use the data to extract 
the ratio C( 2 Vf/v2 using Eq. (4). Using the YBCO6.95 
data, we obtain a velocity anisotropy vf/v2 = 6.8 if we 
assume that a — 1. [|7| 

Alternatively, by comparing the measured slope 
dp s /dT of the YBCO6.95 and YBCOg.60 samples, we 
see that the slopes are almost the same, showing that 
a 2 VF /v2 is almost independent of doping. From tunnel- 
ing data we know that the maximum gap Ao slightly 
increases with underdoping. || This implies that a is al- 
most independent of x. This is the experimental evidence 
that the Fermi liquid scenario a = x does not apply to 
the underdoped cuprates. 

It is useful to compare Eq. (4) with the standard BCS 
expression which is usually written in the form [pl[ 



p s (T) = Ps (0) 1 



(21n2)T 
A 



p s (0)-p s (0)(21n2)T/A 



(8b) 



This expression does not include Fermi liquid correction 
and should be compared with Eq. (6) with a — I. We 
note that in BCS theory p s (0) is independent of x and the 
second term in Eq. (8) is in exact agreement with the sec- 
ond term in Eq. (6) as it should be, because the derivation 
leading to Eq. (6) is completely general. The first terms 
in Eq. (8) and Eq. (6) do not agree because the standard 
BCS theory does not apply to a doped Mott insulator 
and does not include the physics leading to a spectral 
weight proportional to 1. It is clear that this feature of 
Eq. (8) does not agree with experiment on underdoped 
cuprates. If one ignores this and fits the normalized data 
p s (T)/p s (0) to Eq. (8a), one would reach the incorrect 
conclusion that the energy gap Ao is proportional to T c 
in underdoped cuprates. |52| We emphasize that Eq. (6) 
includes Eq. (8) as a special case and must be used in 



place of Eq. (8) for a correct analysis of the data. 

We can also estimate the size of the vortex core using 
this picture. The idea is to identify the core size as the 
point where the critical current is reached. If we replace 
— eA/c in Eq. (2) by the gauge invariant superfluid ve- 
locity v s = i (V6> — ^fA), we see that the quasiparticle 
energy shifts up or down in the presence of v s and quasi- 
particles are generated at the Fermi energy, contributing 
to a normal fluid density. Near the vortex core, v s grows 
as 1/-R, so that the normal fluid density grows and even- 
tually drives the critical current to zero. This allows us 
to estimate the core size to be 



Ri 



Vf 



xirA {t/3J) ttT c 



(9) 



Note the factor x appears in the denominator. We note 
that in BCS theory, the coherence length can be written 
either as vf/kTc or vf/Aq. The two forms are equiva- 
lent because the ratio 2Ao/fcT c is a constant. In our case 
this ratio depends on x and it is not clear a priori which 
form is correct for the coherence length. Equation (9) 
shows that vf/ttT c is the correct form for the coherence, 
and not vf/ttAq. One consequence of this is that the un- 
derdoped cuprates are in fact not short coherence length 
superconductors. Jlq ] The number of holes per coherence 
volume actually grows as x~ l with decreasing doping. A 
second consequence is that H C 2 (due to orbital effects) 
is predicted to scale as x 2 . Within this picture it is also 
clear that in underdoped cuprates the state inside the 
vortex core should retain the large gap Ao, just as the 
normal state above T c . 

We can now estimate the condensation energy using 
the relation AE = H 2 /8tt and H 2 = H cl H c2 . Not- 
ing that H c i is proportional to p s (0)/m w xt while 
H C 2 = <t>o/R 2 is proportional to x 2 , we find that AE 
is proportional to x 3 , i.e., 



AE « x(T 2 /J)(t/3J) 



(10) 



This is in contrast to the BCS expression AE w T^/ep- 
Equation (10) also follows from a picture where only 
the quasiparticles with energy less than T c are affected 
by the transition to the normal state. The area of the 
Brillouin zone occupied by these excitations is of order 
(T c /J)(T c /Aq), so the total energy change per area is 
of order T 3 /JAo which agrees with Eq. (10). Thus even 
when expressed in terms of T c the condensation energy is 
much less than the BCS value in the underdoped system. 
There is evidence for this suppression of the condensation 
energy from specific measurements. H 



IV. THE SU(2) FORMULATION OF THE T-J 
MODEL 

We now return to discuss the microscopic theory. 
While the mean field phase diagram is in qualitative 
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agreement with experiments, the £7(1) formulation suf- 
fers from a number of deficiences if we try to improve 
the mean field theory by including gauge fluctuations at 
the Gaussian level. In the spin gap phase the problem 
lies with the fact that the MF theory is a pairing the- 
ory of fermions and carriers with it some features of su- 
perconductivity. For example, the gauge field is gapped 
by the fermion pairing via the Anderson-Higgs mecha- 
nism. This leads to a reduction of gauge fluctuations 
which actually destabilize the pairing phase. |33) A sec- 
ond problem is that if we introduce residual interaction 
between the fermions and bosons to form an electron, 
the electron spectrum will always have nodes. This is 
because the node structure in the pairing state is tied to 
the Fermi level and is very resilient to interactions. Thus 
we have difficulty reproducing the "Fermi surface seg- 
ments" which are apparently observed in photoemission 
experiments. In the superconducting phase we have con- 
densation of the bosons and the quasiparticles become 
well defined. While this feature is in agreement with ex- 
periment, the current carried by the quasiparticles turns 
out to be reduced so that in Eq. (2), a = x. As we have 
seen, this leads to a serious disagreement with the doping 
dependence of the temperature coefficient of the London 
penetration depth. In order to circumvent these difficul- 
ties, we were led to a new formulation of the t-J model 
which is designed to be more accurate near half filling. 
We briefly outline the SU(2) formulation below. HQ 
In this new formulation we introduce an SU{2) doublet 
of boson fields b T — (61,62)) in addition to the fermion 
doublet ipi = (i/'f, ip\)- The physical electron is repre- 
sented by the SU(2) singlet formed out of these two dou- 
blets, c-f = -^b T ip, Cj = where ip — ir 2 tp*. We 
are motivated by the observation made by Affleck et al. 
p6| that at half-filling (x — 0) the fermion representa- 
tion of the t-J model has the SU(2) symmetry in that a 
spin-up electron can be represented by a spin-up fermion 
or the absence of a spin-down fermion. In the £7(1) for- 
mulation this symmetry is broken as soon as x 7^ 0, and 
out of a infinte degeneracy of states, the ci-wave fermion 
pairing state is picked out as the MF solution. In con- 
trast, even at the mean field level, the low lying states 
which are missing in the £7(1) mean field theory are in- 
cluded in the new SU(2) formulation. For example, the 
spin gap state can be described equally well as the d- 
wave pairing state, or a staggered flux phase, where the 
fermions see gauge fluxes which alternate from plaquette 
to plaquette. The SU(2) gauge transformation relates 
these states and guarantees that there is no breaking of 
the translational symmetry. The fermion spectrum ex- 
hibits a d-wave type gap, with maximum gap at (71", 0) 
and nodes at (7r/2,7r/2). We compute the physical elec- 
tron spectral function, which at the mean field level, is 
a convolution between the fermion and boson spectra. 
We further introduced a residual interaction between the 



fermions and bosons. The resulting spectra can be com- 
pared with photoemission experiments and have the fol- 
lowing features. The spectra consist of a coherent part 
with spectral weight x and dispersion of order J and 
a broad incoherent part. The coherent part closely re- 
sembles the fermion dispersion. The residual interaction 
broadens and shifts the nodes at (ir/2,ir/2) so that we 
obtain a "Fermi surface segment" near (tt/2,tt/2). Away 
from this segment a gap appears in the excitation spec- 
trum which grows to its maximal magnitude near (0, n). 
This behavior is in qualitative agreement with the angle- 
resolved photoemission experiment. ||,[| 

We have also studied the fermion spectrum and how it 
is affected by gauge fluctuations. We found a logarithmic 
correction to the fermion velocity and we successfully fit- 
ted the magnetic susceptibility and the specific heat in 
the spin gap state. |37j 

In the superconducting state we need to address the 
issue of the current carried by the quasiparticles. To ex- 
pand on this point further, we note that in the original 
£7(1) gauge field formulation of the t-J model, the predic- 
tion for p s (T) takes the form of Eq. (6) with a — x and 
therefore is in strong disagreement with experiment. This 
follows simply from the Ioffe-Larkin rule which states 
that the inverse of the response function of the fermion 
and boson should add to give the physical inverse re- 
sponse. In the superconducting state, the fermion and 
boson acquire superfluid densities pf and p s so that 

p- 1 {T)=p F \T)+p B l {T) (11) 

where pf ~ (1 — x) and pb ~ x. However, only the 
temperature dependence of pp depends on the qp gap 
structure and is expected to be of the form pf(T) ~ 
(1 — x)(l — T/Ao), whereas the temperature dependence 
of ps arises only through the excitation of sound mode 
and should be higher power in T, which can be ignored. 
Inserting these into Eq. (11) we see that p s (T) is pre- 
dicted to be x— x 2 T/Aq. Basically in the £7 (1) gauge the- 
ory the mismatch of the Fermi surface area and the Drude 
spectral weight (or p s in the superconducting state) is 
solved by a Landau parameter, so that a = x. Thus we 
may conclude that it is not sufficient to treat the gauge 
fluctuation only to quadratic order as in the Ioffe-Larkin 
theory. 

We believe this difficulty is tied to the notion of Bose 
condensation as a way of achieving superconductivity. 
The reason is the following. The electron operator Ck is a 
convolution of the fermion and boson operator in momen- 
tum space. Let us suppose that the external A field cou- 
ples only to the boson (this is true in the SU (2) formu- 
lation and is approximately true in some gauge choice in 
the £7(1) formulation). In the presence of A, 6 q — > 6 q +A 
so that after the convolution Ck — > Ck+A and ek — > £k+A 
as expected. Thus jk = —ede/dA = —ede/dk. Let us 
see what happens in the superconducting state. If wc 
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assume that the fermions are already paired, supercon- 
ductivity can be driven by the condensation of bosons 

< &k=o >= 0. However, in the presence of A, the Bose 
condensate remains rigid and stays in the k = state. 
This is clearly seen in the Ginsburg Landau theory for 
the free energy | (V — 2eA/c) b\ 2 where < 6k=o >7^ in 
the presence of A is responsible for the Higgs mechanism 
and the London penetration depth. Upon convolution, 
we see that for the electron operator, k is not shifted by 
A so that e(k) is independent of A. The qp now carries 
no current! In the U(l) formulation, the gauge field a 
causes a small shift in the Fermion spectrum and leads 
to Eq. (2) with a — x. This is clearly an unacceptable 
situation and can be seen most acutely for the qp at the 
Fermi surface along the (7r,7r) direction. Here the energy 
gap vanishes so that the qp in the superconducting state 
is basically the same state above T c . Yet, according to 
the Bose condensation scenario, the current carried by 
this qp drops abruptly below T c . 

Now that we have identified the problem, we can see 
that there are two possible ways to avoid it. The first is 
to argue that due to fluctuations, only a small fraction 
of the bosons are in the condensate and we can reduce 
the problem, but not eliminate it. We call this the single 
boson condensation (SBC) scenario. The result is that 
a can lie anywhere between x and 1, and most likely 
somewhere in between. A second possibility is allowed in 
the SU(2) formulatin but not in the U(l) formulation. 
In SU(2) theory there are two species of bosons b\ and 
62 and we can pair them to form a gauge singlet pair 

< 61 (i)&2 (j) >7^ 0. We shall call this the boson pair con- 
densation (BPC) scenario. Since < b\ >—< hi >= 0, the 
problem is avoided and we find that a = 1 . This is really 
a consequence of continuity because in this scenario the 
superconducting qp along is smoothly connected to 
the electron state above T c . This result comes out of an 
explicit calculation which we outline below. 

In SU(2) theory we go beyond MF theory by calcu- 
lating the electron propagator through a ladder diagram 
fj4| , j35f to include effects of pairing between the boson 
and the fermion. Here we will consider only the simplest 

on-site interaction V ^c|c| + c| c i) > which? when written 
in terms of bosons and fermions, generates an attraction 
between bosons and the fermions if V > 0. There are also 
other pairing interactions, but they will not modify our 
results qualitatively. The resulting electron propagator 
is given by 



G A (w,k) 



-z(c T cj) 
-i(c\c\) 



' l \ c t c l) 
-«(4 c l) 



G .a(w, k) 



*b,A(w,k) 
-Go,a(-w, -k) 



(12) 
(13) 



We first consider the second scenario where there are 
no SBC, but there is a nonzero Fq.a proportional to 
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FIG. 4. Schematic illustration of the QP dispersion (the 
pole location of G for (a) normal state, and SC state with (b) 
< x pc < x and (c) x pc = x. The line thickness indicates 
the size of the residue of Gu, and the dashed line indicates 
vanishing residue. The momentum scan is along the straight 
line in (d), where the curved segment is the FS segment in 
the normal state. 



the boson pair parameter x pc . For A = 0, the poles 
of G ! n(o;,k) comes in pairs of opposite signs, just as in 
BCS theory. However the total residue is 2 (i-vg 
significantly reduced from the BCS value. There are two 
positive branches which determine the qp excitations 



E^ c \k) 



where 



-E± = ±l/(e-(i) 2 



/' 



(14) 



(15) 



and fi = ~ 4(i_fy-G. ) ■ I n order to interpret those results, 
let us first consider the normal state which is recovered 
by setting x pc = in Eq. (14) and Eq. (15), yielding 
the normal state dispersion E± = E±(x pc = 0). This 
corresponds to a massless Dirac cone initially centered 
at (±7r/2,±7r/2) when V = which is the MF fermion 
spectrum of the staggered-Flux (s-Flux) phase. The ef- 
fect of V (the boson-fermion pairing) is two-fold. The 
/2 inside the square-root shift the location of the node 
towards (0,0) by a distance Afc = — jl/vp while the last 
term shift the spectrum upwards. The cone intersects the 
Fermi energy to form a small Fermi pocket with linear 
dimension of order x. As shown in Fig. 4(a), the spectral 
weight is concentrated on one side of the cone, so that 
only a segment of FS on the side close to the origin car- 
ries substantial weight. This is the origin of the notion 
of "FS segment" introduced in Ref. 

Now let us see what happens in the SC state when 
Xpc 7^ 0. Equation (14) takes the standard BCS form if 
E± is interpreted as the normal state dispersion. How- 
ever, E± differs from the normal state spectrum E± by 
the appearance of the term — (x pc A/x) 2 in Eq. (15). 
Close to the node this term is small so that qualitatively 
the spectrum develops from the normal state in a BCS 
fashion, as shown in Fig. 4(b). This is particularly true 
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if the higher energy gap between the two branches is 
smeared by lifetime effects. Thus we see that the "FS 
segment" is gapped in a BCS-like fashion. However, the 
velocity vi in the (1, —1) direction, being proportional to 
Xp C /x, does not extrapolate to the gap at (0, 7r) (which is 
essentially independent of x pc ), but crosses over to it at 
the edge of the FS segment. It is worth remarking that 
in the special case x pc — x, reduces to the standard 

BCS form with the normal state dispersion e(k), a chem- 
ical potential 2jl and a SC gap A(k). The high energy 
gap closes and spectral weight on one branch vanishes, 
yielding a BCS spectrum as shown in Fig. 4(c). 

We have also calculated the effect of constant A on 
the qp dispersion, to linear order of A. This adds a term 
-j± • A to Eq. (14) where j± is interpreted as the cur- 
rent carried by the qp. We recall that in standard BCS 
theory, the current is given in term of the normal state 
spectrum by c<9a£a = e<9k£ because £A(k) = e(k + |A). 
Remarkably this is almost true in our case in the sense 
that j± is given by c8a.E±^ where E± t A is obtained by 
replacing k by k+ |A in e, jl and A everywhere in Eq.|l5| 

except for the term (^A) , which is kept independent 
of A. Near the node, A is negligible so that the current 
is very close to e&^E ~ ed-^E N (which becomes exactly 
e<9k£ along the diagonal), thus reproducing Eq. [|. We 
have checked numerically that even away from the node 
in the region of the "FS segment" , the current is remark- 
ably close to edkE N , which can be quite different from 
the BCS value ed^e near the edge of the FS segment. 

From Eq. (6), the temperature dependence of the 
London penetration depth gives a direct measurement 
of a 2 -^ 1 - Density of states measurements using the T 2 
coefficient of the specific heat yields vfV2- The Fermi 
velocity can be estimated from transport measurements 
or high resolution photoemission experiment. Thus in 
principle the quantities a, vf and vi can be measured. 
It is of course of great interest to establish how close a 
is to 1, or whether v-i is reduced with respect to that 
extrapolated from the energy gap at (0, 7r) measured by 
photoemission or tunneling. Crude estimates made in 
Rcf. J27j suggest that a is consistent with 1 but a more 
precise measurement is clearly called for. 

Finally we comment on finite temperature behaviors. 
In addition to the reduction of superfluid density due 
to thermal excitation of qp, we expect x pc to decrease 
with increasing T, leading to a reduction of V2- V2(T) = 
a PC (o) ^ 2 (0)- ^ s T reaches T c , x pc = v 2 = and the nodes 
of i?( sc ) become the "FS segment" while the spin gap near 
(0, 7r) remain finite. We see that x pc plays the role of the 
order parameter of the transition, so that we may expect 
the temperature dependence of x pc to be described by 
a Ginzburg-Landau theory with X-Y symmetry near the 
transition. 



V. CONCLUSIONS AND OPEN ISSUES 

We believe the SU(2) slave boson theory captures the 
basic physics of the underdoped cuprates. The many 
anomalous properties associated with the spin gap forma- 
tion are explained in a natural way. Superconductivity 
with d-wave pairing symmetry emerges naturally, with 
quasiparticle excitations which are remarkabley similar 
to BCS theory. However, the microscopic mechanism is 
completely different in that the SC state is not formed 
out of pairing of normal state quasiparticles via exchange 
of some effective interaction. Instead, it is the coherence 
of the charge degrees of freedom which converts the spin 
gap phase to the SC state. Many open issues remain, 
however, and we list a few of them below. 

1) Our discussion of the electron spectrum in the nor- 
mal state is still at a crude level. We treat the 
bosons as "nearly" bose condensed with a relatively 
narrow spectral function. Thus we do not have a 
theory of the lineshape. One of the most important 
features of the photoemission experiment is that a 
narrow qp peak forms out of a broad lineshape as 
the SC state develops out of the normal state. We 
are unable to describe this evolution at present. 
A narrow spectral line is very natural in the sin- 
gle boson condensation scenario but not as obvious 
in the boson pair condensation scenario. Thus we 
have not achieved a quantitative description of the 
recombination of spin and charge to form quasipar- 
ticles in the superconducting state. A related issue 
is that in our theory the spin gap state and the 
SC state share the same energy scale, i.e., the en- 
ergy gap A at (0,7r). Empirically A J/3, in 
rough agreement with the gap calculated in mean 
field theory. Recently, Shen and collaborators |||] 
have focused on a higher energy scale (of order J 
to 2 J) which characterizes the location of the peak 
in the ARPES spectrum, and argued that it is the 
peak energy which is smoothly connected with the 
insulator at half-filling. In this scenario one would 
need a separate mechanism to produce the leading 
edge shift and the SC energy gap. In our scenario 
we have only one energy scale Ao and the burden 
upon us is to show that the lineshape may exhibit 
a peak at high energy of order J. 

2) We do not have a satisfactory theory of the trans- 
port of the normal state. This is related to the still 
lack of understanding of how the spin-charge sepa- 
ration state in the normal state evolves to the well 
defined qp in the SC. We can only provide a phe- 
nomenological picture of gradual binding between 
holons and spinons to form physical holes as the 
temperature is decreased. 
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3) The mean field theory underestimates the spin fluc- 
tuation near (ir, w). While inclusion of gauge fluc- 
tuations leads to a satisfactory fit of the specific 
heat and uniform spin susceptibility, p3]| it is ex- 
pected Jill that gauge fluctuations will strongly en- 
hance the spin fluctuation near (tt, tt) but detailed 
calculations have not been carried out. This strong 
enhancement is needed to explain the strong peak 
in the Cu NMR relaxation at a temperature T* 
which is low compared with the spin gap energy Aq. 
As an intermediate step, we recently carried out a 
RPA calculation of the spin fluctuation near (tt, it). 
|Q By tuning a single parameter (the effective ex- 
change coupling in RPA) we are able to account 
for the resonance peak seen in neutron scattering 
in the SC state and its evolution with reduced dop- 
ing. [ fl3"| However, at present we cannot explain the 
neutron scattering and the copper NMR within the 
same RPA theory. 

The work reviewed in this paper has been done in close 
collaboration with X.G. Wen and I have benefitted from 
collaboration over the years with N. Nagaosa, T.K. Ng, 
Derek K.K. Lee, Don H. Kim, and J. Brinckmann. This 
work was supported by NSF through the MRSEC pro- 
gram DMR 98-08941. 
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